A formalism for the construction of binary neutron stars with arbitrary circulation 
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Most numerical models of binary stars - in particular neutron stars in compact binaries - assume 
tiie companions to be either corotational or irrotational. Either one of these assumptions leads to a 
significant simplification in the hydrodynamic equations of stationary equilibrium. In this paper we 
develop a new formalism for the construction of binary stars with circulation intermediate between 
corotational and irrotational. Generalizing the equations for irrotational fiow we cast the Euler 
equation, which is an algebraic equation in the case of corotational or irrotational fiuid fiow, as an 
elliptic equation for a new auxiliary quantity. We also suggest a parameterized decomposition of 
the fiuid fiow that allows for a variation of the stellar circulation. 



PACS numbers: 



I. INTRODUCTION 



Significant effort has gone into the numerical modeling 
of compact binaries containing neutron stars (see, e.g., 
[ij for a review). Quasiequilibrium models of such bina- 
ries usually assume the existence of an approximate he- 
lical Killing vector that generates the binary orbit. The 
existence of such an approximate Killing vector implies 
that, in a coordinate system corotating with the binary, 
the binary companions appear approximately stationary. 
However, this notion leaves the rotation, or circulation, 
of the individual stars undetermined. For example, each 
star could always show the same side to its companion, 
making it corotational, or it could have zero circulation 
as seen in an inertial frame, making it irrotational. In 
principle, of course, each star could be rotating arbitrar- 
ily, about an arbitrary axis, and potentially even differ- 
entially. 

As long as neutron stars can be modeled as ideal flu- 
ids, constructing neutron star solutions amounts to find- 
ing solutions to the Euler equation and the continuity 
equation. For equilibrium solutions, the Euler equation 
simplifies in two special cases, namely for corotational 
and for irrotational binaries. In both cases the Euler 
equation reduces to a total derivative and can be inte- 
grated once to yield an algebraic equation for the fluid 
enthalpy. 

For corotational binaries the continuity equation is 
satisfied identically, which simplifies the resulting equa- 
tions significantly. Not surprisingly, this assumption was 
adopted in the first quasiequilibrium studies of both bi- 
nary neutron stars (see, e^g., 0, Q) and black hole- 
neutron star binaries (e.g. [^). Corotational flow does 
not seem very realistic, however, since tidally locking a 
neutron star in a compact binary would require an un- 
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physically large viscosity (see [1, @|). 

Instead, it seems more realistic to assume the fluid to 
be irrotational, which is a reasonable assumption as long 
as the orbital frequency is much larger than the spin fre- 
quency that the neutron stars would have in isolation. In 
this case the velocity may be written as the gradient of a 
velocity potential. The continuity equation then becomes 
an elliptic equation for this velocity potential. This for- 
malism was developed independently by 0, H, H, 
(see also [llj for a comparison between the different ap- 
proaches) , and numerical solutions have been constructed 
both for binary neutron stars (e.g. [13, 1X3, 14], sec also 
for Newtonian models) and black hole-neutron star 
binaries (e.g. EESIIl). 

Constructing binaries containing rapidly rotating 
stars, e.g. millisecond pulsars, requires a formalism that 
allows for arbitrary spin configurations. In this case the 
Euler equation no longer yields a trivial first integral. So 
far, neutron stars with arbitrary circulation have been 
constructed only by [ll|, hereafter MS. MS make several 
assumptions in their treatment of neutron stars with ar- 
bitrary circulation. Instead of using the Euler equation, 
they solve the Bernoulli equation, which states that a 
certain quantity is constant along fluid trajectories. MS 
then find approximate solutions by assuming that this 
quantity is constant globally, which turns out to be cor- 
rect for corotational and irrotational fluid flow, but not 
in general. While this approximation may be quite ac- 
curate, it is of interest, at the very least as a point of 
principle, to identify a formalism that avoids this approx- 
imation. MS further assume that a solenoidal part of the 
fluid flow is centered on the stellar centers, determined 
by the point of maximum baryonic density. This may 
also be a good approximation, and yields corotational 
and irrotational models in good agreement with other 
approaches, but again it would be desirable to develop a 
formalism that avoids this assumption. 

In this paper we develop a formalism that improves 
on that of MS in several ways. We suggest a way in 
which the Euler equation can be solved for circulation 
intermediate between corotational and irrotational, and 
specify the fluid flow in a geometrically motivated fashion 
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that does not make any assumptions about the centers of 
rotation. After stating some prehminaries in Section |TT] 
we will, for the sake of clarity, first develop this formalism 
in a Newtonian framework in Section IIIII We will then 
treat the problem in a relativistic framework in Section 
IIV( and will briefly summarize in Section|Vl In Appendix 
|X]we will briefly discuss an alternative decomposition of 
the fluid flow. 



II. PRELIMINARIES 

We adopt a 3+1 decomposition of the spacetime metric 
and write the spacetime metric gat in the form 

gabdx'^dx'' = -a^dt^ + -f^jidx' + P'dt){dx^ + f3^dt). (1) 

Here a is the lapse function, is the shift vector and 
7ij is the spatial metric induced on a spatial slice S of 
constant coordinate time t. The normal on this slice is 

n- = a-\l,-(3'). (2) 

Here and in the following we denote spacetime indices 
with letters a,b,c, . . . and spatial indices with letters 
fc, . . .. The coordinate congruence t"", which connects 
points with the same spatial coordinate labels on differ- 
ent slices S, is then given by = an" + 

We also assume the existence of a helical Killing vector 
^{Jg[, and assume it to be normalized in such a way that 
its time component is unity, = 1. We can write this 
Killing vector as 

^Li = i° + k\ (3) 

where the vector is purely spatial and describes the 
orbital rotation 

In cartesian coordinates, for a rotation about the 2:-axis, 
we may also write the spatial components of A;" as 

^n{-y,x,0). (5) 



Here is the fluid velocity, po the rest-mass density, P 
the pressure, and (/)n the Newtonian potential, which has 
to satisfy the Poisson equation 

D^^N = 4^Po. (8) 

We denote a covariant derivative with D^, and the 
Laplace operator with D'^ = D^Di. It is also useful to 
define as the fluid velocity relative to fc°, 

v' = k' + V\ (9) 

i.e. y is the velocity as seen in a frame corotating with 
the binary. For isentropic flow we also have 

p^^D,P ^ D,h. (10) 

Here h is the specific enthalpy, which, adopting a conven- 
tion more common in a relativistic framework, is defined 
as 

/7^= 1 + e + P/po, (11) 

where e is the specific internal energy density. We also 
assume that we are given an equation of state that relates 
the enthalpy h to the density po- 

For stationary solutions the Lie derivative of the hy- 
drodynamic variables po a-nd has to vanish, i.e. 

'^Choi^i = dtVi + C\^Vi = dtVi + kWjVi-VjDik^ = (12) 

and similar for pQ. Combining equation (|12|) with the 
Euler equation ^ we then obtain 

A (^h - + + + V^DjV, - D,vj) = 0, 

(13) 

while the equation of continuity becomes 

D^{poV') - 0. (14) 

Equations (fT5|) and form the Newtonian equations 
of stationary equilibrium. 

B. Corotational and irrotational flow 



III. NEWTONIAN TREATMENT 

A. Newtonian equations of stationary equilibrium 

The Newtonian equations of hydrodynamics are the 
Euler equation 

dtVi + vWjVi = -pQ^DiP - Di(j)^ (6) 

and the continuity equation 

5tpo + A(po«') = 0. (7) 



It is easy to see how the above equations simplify for 
corotational or irrotational fluid flow. 

For corotational flow, we must have ~ 0, in which 
case the continuity equation (|14p is satisfied identically. 
Furthermore, the second term in the Euler equation (jl3p 
vanishes identically, so that the first term can be inte- 
grated to yield an algebraic equation for the enthalpy ft,, 

- ' + = C. (15) 

Here C is some constant of integration, and can be 
found from Q, which now implies = /c*. Here and 
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below we assume that both stars are identical for sim- 
plicity; the method is easily generalized for non-identical 
stars. 

For irrotational flow, we require the fluid's vorticity to 
vanish, which is automatically the case if we write the 
fluid velocity as a gradient of a velocity potential $, 

= (16) 

In this case, the second term in the Euler equation 
also vanishes identically, so that we again obtain an al- 
gebraic equation for the enthalpy h, 

h-^v^ +VjV^ +(bN = C. (17) 

The continuity equation (|14p now becomes an elliptic 
equation for the velocity potential <i>, 

1)2$ = (fc* -D*$)Aln/9o- (18) 

To ensure that the right hand side remains finite at the 
stellar surface, where po becomes zero, we must have 

(fc'-n'$)Apo|_f = 0. (19) 

Given that DiPQ is normal to the stellar surface, this 
relation furnishes a Neumann boundary condition for $. 
We could have derived the same condition by requiring 
that, at the stellar surface, the fluid flow be tangent to 
the surface. 



C. Intermediate circulation 

For intermediate circulation, the second term in the 
Euler equation no longer vanishes, and we can no 
longer find a simple first integral. One approach to find 
a solution is to take a divergence of the equation, which 
yields an elliptic equation 

D'^H + {V\DjV, ~ DiVj)) = (20) 

for a new auxiliary quantity 

H = h-^v'^ + VjV^ +(1)^. (21) 

In the following we will interpret this as an equation for 
the specific enthalpy h, 

h = H +^v^ -VjV^ - (1)^. (22) 

Before proceeding we need to specify the circulation of 
the fluid. To do so, we choose to decompose the velocity 
Vi according to 

V, ^ D,<S> + 7^k„ (23) 
which, from ([9]), implies 

V' = D'<^> + {f]~ l)k\ (24) 



Here $ is a generalized velocity potential, and rj is an 
arbitrary constant. As we will discuss in Section IIII Dl 
we recover irrotational flow for ry = and corotational 
flow for 77 = 1. With our decomposition t] there- 

fore parameterizes a sequence of stellar models connect- 
ing corotational and irrotational models. By generalizing 
the decomposition it should be possible to change 
the axis of the internal rotation, or change the degree of 
differential rotation. 

Inserting the decomposition into equation ([^0]) we 
first note that the velocity potential $ again drops out 
of the term 

DjVt ~ DiVj = T] (Djki - Dikj) = r]Kij, (25) 

where we have defined for convenience. Evidently, Ky- 
is antisymmetric by construction and can be computed 
directly from ([5]) . The second term in equation (|20p then 
becomes 

(26) 

Inserting to evaluate the first term on the right hand 
side we find 

(DV^)Ky = iDWi^+ir]^l)D'k^)K,j 

= {'n-i){D^y)K,„ (27) 

where the second equality holds because D^D^^ is sym- 
metric while Kij is antisymmetric. That means that all 
second derivatives of $ disappear in equation (|20p . leav- 
ing us with 

D'^H+Tj{Tj-l){D'k^)Kij+rj {D^^ + {v- ^W) D'nij = 0. 

(28) 

This is the form of the Euler equation that generalizes 
equations (jlSp or (|17p for intermediate circulation. In- 
stead of an algebraic equation for the enthalpy h we now 
have to solve an elliptic equation for the new auxiliary 
quantity H . Given solutions for the other fields we can 
then find the enthalpy h from ((221) • On the stellar surface 
the density vanishes, meaning that h approaches unity. 
This requirement leads to the Dirichlet boundary condi- 
tion 

Hl^,,^l-\v^+v,V= (29) 

on the stellar surface. 

To find an equation for the velocity potential $ we 
insert ((24|) into the continuity equation (fTi]) . which yields 

1)2$ ^ - (1)'$ + {n- 1)F) A In Pq. (30) 

To ensure that the right hand side remains finite on the 
stellar surface we must have 

(i^'$ + (77-l)/c^) ApoL,,f = 0. (31) 

As in the case of irrotational flow, this conditions forms 
a Neumann boundary condition. 
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We have now assembled all equations and boundary 
conditions needed to construct binary stars with circula- 
tion intermediate between corotational and irrotational 
flow in a Newtonian framework. In addition to the Pois- 
son equation ([8]), we need to solve equation ((28|) . subject 
to the boundary condition (I29p . for iJ, and equation ([30|) . 
subject to boundary condition (|3ip. for Given these 
solutions we can find Vi and from equations (|23p and 
(|24)) . and h from ((22)) . Finally, po can be computed from 
/i given an equation of state, and O is determined by 
enforcing a circular orbit (see, e.g., fisj). 

D. Limiting cases 

Before proceeding to a relativistic treatment of the 
problem it is useful to see how the equations for coro- 
tational and irrotational flow appear as limiting cases of 
the equations for intermediate circulation. Clearly, our 
equations for intermediate circulation are consistent with 
equations and (HH), and are therefore satisfied both 
for corotational flow with ~ and irrotational flow 
with Vi = Di<^. In the following we demonstrate that 
these solutions are also consistent with the decomposi- 
tion ((231) . 

We start with irrotational flow, which we obtain for 
?7 = 0. In this limit, the decomposition (|23p evidently 
reduces to ((TC|) , and both the equation ((50)) and boundary 
condition ()3ip reduce to their equivalents for irrotational 
flow, equations ((T8)) and ((T9)) . For 77 = 0, equation ((28)) 
reduces to a Laplace equation for iJ, which is uniquely 
solved by iJ = Ch, where Ch is a constant, when H = 
Ch on the boundary. This is consistent with the other 
equations, which can be seen by evaluating equation ITT)) 
on the stellar surface. As a result, the equation for the 
enthalpy h reduces to equation ((T7)) . as expected. 

We obtain corotational flow by setting 77 = 1. To see 
this, we first observe that both the equation ((30)) for the 
velocity potential $ and its boundary condition ((^ are 
solved by <& = C$, where C$ is a constant. Equation (IM)) 
then implies 1^* = 0, as expected for corotational flow. 
Equation ((28)) again reduces to a Laplace equation which, 
as in the irrotational case, is solved by i? = Ch- The 
equation for the enthalpy h therefore reduces to equation 

IV. RELATIVISTIC TREATMENT 

A. Relativistic equations of stationary equilibrium 

The relativistic equations of hydrodynamics follow 
from the conservation of energy-momentum, 

VbT"^ = 0, (32) 

whose spatial components form the relativistic Euler 
equation, and the conservation of rest-mass 

Va(pow") = 0. (33) 



Here T°'^ is the fluid's energy-momentum tensor, its 
four-velocity, and Va denotes a covariant derivative asso- 
ciated with the four- metric gab- Following the approach 
of [l^l we can write u° as 

= + (34) 

where the helical Killing vector S^^^^ was introduced in 
equation ([3]) and where V"" is purely spatial, UaV^ — 0. 
Just like its Newtonian counterpart in equation ([5]), y° 
now measures the fluid velocity in a frame comoving with 
the binary. It is also convenient to define 

Ui = ji^hua- (35) 

In terms of the normalization of the four-velocity 
Uqu" = — 1 may be written as 

au' = {1 + h-^ u,u'y^\ (36) 

For stationary solutions the Lie derivative of the hy- 
drodynamical variables along ^jjgj again has to vanish. 
As shown in [Toj , the relativistic Euler equation then be- 
comes 

A (^^ + UjV'^ + V {D,u, - D,u,) = 0, (37) 

while the continuity equation becomes 

A: {au^poV) = 0. (38) 

Here Di is the covariant derivative associated with the 
spatial metric 7^ . 

We also assume that, on the slice E, the spatial metric 
7ij, the lapse a, and the shift /3* are constructed by solv- 
ing some subset of Einstein's equations, for example by 
using the conformal thin-sandwich decomposition. These 
equations then play the role of the Poisson equation ([8]) 
in the Newtonian framework. As before we further as- 
sume that we are given an equation of state that relates 
the specific enthalpy h to the rest density po- 



B. Corotational and irrotational flow 

As in the Newtonian case it is easy to see how the 
equations simplify for both corotational and irrotational 
flow. 

For corotational flow the four-velocity u° is aligned 
with the helical Killing vector ^^^j, so that y° = 0. The 
only remaining term in the Euler equation (|37p can then 
be integrated to yield 

-^=C, (39) 

where C is some constant, while the continuity equation 
((55)) is satisfied identically. 
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To obtain irrotational flow, with vanishing relativistic 
fluid vorticity, we now assume 



(40) 



where $ is a velocity potential (see [ifll for a discussion) . 
With this choice the second term in the Euler equation 
P7p vanishes again, and we can find the first integral 



+ u^V' = C. 



(41) 



Combining equations ((34)l . ([36)1 and (|4T|) we also find 
1 



ah 



(42) 



where we have defined a rotational shift vector 

B'' = 13'' + k" . (43) 
Using ^ and — an"- + j3°- we may also write 5° as 



an 



(44) 



We will derive a more general version of equation (f42|) 
for intermediate circulation below. 

Note that can always be expressed as 



(45) 



For irrotational flow we can now insert (I40p . 



y = ^-D'$ - B\ 

and insert this result into the continuity equation 
obtain 

D,{apah-^D'<P) - Di{au^poB') = 0. 

We may eliminate u* from this equation by inserting 
which yields 



DH^D,i^±^B- 



(46) 
to 

(47) 



^±^^B^-D^^]Dan^. (48) 



As expected, this equation reduces to in the Newto- 
nian limit. A Neumann boundary condition again results 
from requiring that the right hand side remain finite at 
the stellar surface. 



^±^B^-D^^]Dr' 



= 0. 



(49) 



To complete the argument, we could insert equations 
(liO)) . (li^ and (gni) into (gT]) to obtain an equation for 
h in terms of and independently of w'. 



C. Intermediate circulation 

As in the Newtonian case, we observe that for interme- 
diate circulation the second term in the Euler equation 
([57)1 no longer vanishes, so that we can no longer find a 
first integral. Taking a divergence of the equation, we 
can again obtain an elliptic equation 

D^H + D' {V^iDjUi - Diiij)) = (50) 

for a new auxiliary function 



H=—+ u,V' 



(51) 



We will again interpret this as an algebraic equation for 
the specific enthalpy h. 

Generalizing the Newtonian expression (|23p we decom- 
pose Ui according to 

u,^D,^ + r]ki. (52) 

Repeating the steps in equation (I46p we now find 

1 



-{D'^ + Tjk')~ B' 



(53) 



We point out that there are many different ways of gen- 
eralizing the Newtonian decomposition (j23p . and we will 
discuss an attractive alternative to the decomposition 
([5^ in Appendix [A1 Moreover, the Newtonian decompo- 
sition represents only one particular sequence that 
connects corotational and irrotational fluid flow; it should 
also be possible to generate models that rotate about an 
arbitrary axis of rotation. Similarly, the relativistic de- 
composition (j52p represents only one possible sequence. 

Before proceeding we derive an equation for m*. Start- 
ing with the deflnition (|5ip we insert equation (|45p to 
obtain 



1 



-it' - B' 



— (1 + h-^u.ii') = H + u,B' 



(54) 



(55) 



We now insert the normalization as well as the de- 
composition ((52)) to find 



au 



{H + B'{i]k,+D,^)) 



(56) 



For the irrotational case, when jy = and H — C, we 
recover equation ([^^ . as expected. 

Returning to the Euler equation (|50|) , we observe that 



(57) 



in complete analogy to the Newtonian case. Expanding 
the terms in (1501) we now find 



1 



D^H + ri{(D^^ + nk^)D'— + -^D'k^ - D'B^ n 



V\-rriD^'^ + ri^) - B' D'nu = 0. 



(58) 
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If desired, equation ([56| can be inserted to eliminate the 
terms Requiring that the density vanish on the stel- 
lar surface, so that h = 1, again leads to a Dirichlet 
boundary condition 



(59) 



We now turn to the continuity equation (j38p . Inserting 
the decomposition ([55)) . we obtain 

A (^(^'* + Vk')) - D,{au'p„B') ^ 0, (60) 
or, eliminating with the help of equation ([51 



A$ - A 



= -i^D^k' (61) 



-A - D'^-rjk' A In 



As before we require that the right hand side remain 
finite at the stellar surface, which yields the boundary 
condition 



surf 



H + B^{rjk,+D,^)^^ ...^n«Po _Q 

(62) 

We have now assembled a complete set of equations 
and boundary conditions that can be used to construct 
relativistic binary stars with intermediate circulation. 
Given a choice for the parameter i] we can solve equation 
([58|l for H and equation (|6ip for $. We can then find the 
combination u^h from (|56|1 . and the velocity expressions 
y and iii from equations ((5^ and ([55]) . Finally, /i can 
be found from the definition (j51[) . po from the equation 
of state, and fl can be determined by enforcing a circular 
orbit (see [l| and references therein). The above equa- 
tions have to be solved together with gravitational field 
equations of the lapse a, the shift /?* (which, together 
with the helical Killing vector determines the rota- 
tional shift vector B^) and the spatial metric 7^^; this can 
be accomplished, for example, by solving the constraint 
equations of Einstein's field equations in the conformal 
thin-sandwich decomposition. All the above equations 
are coupled in complicated ways, and presumably the 
easiest way of constructing numerical solutions involves 
an iterative algorithm. 



D. Limiting cases 

Our equations are based on the Euler equation (|37|) 
and the continuity equation (|38p . which are clearly sat- 
isfied for both corotational flow with ~ and irro- 
tational flow with Ui — Di^. However, it is less clear 
whether these solutions are consistent with our decom- 
position dnm and dnnn- 

As in the Newtonian case we can recover irrotational 
flow for 77 = 0. In this limit, all equations, including a 



constant solution H = Ch for H , are consistent with the 
corresponding equations for irrotational flow in Section 
lIVBl 

On the other hand, we do not seem to recover exactly 
corotational flow with y = for 77 = 1 for the decom- 
position (jSSp . since it is not clear how, for 77 = 1, we 
would obtain a solution for $ that exactly cancels the 
other terms in ([55]) . As we have seen in Section UlIBi we 
do obtain corotational flow for 77 = 1 in the Newtonian 
limit; for relativistic solutions, deviations from corotation 
can therefore be at most at the post-Newtonian order. 
Our decomposition (j52p therefore represents a sequence 
that connects irrotational flow with approximately coro- 
tational flow. In Appendix |^ we present an alternative 
decomposition of the velocity field, that does reproduce 
exactly corotational fiow for 77 = 1. However, this de- 
composition also introduces several additional derivative 
terms that may prove problematic, as we discuss in Ap- 
pendix \K\ 



SUMMARY AND DISCUSSION 



In this paper we propose a formalism for the construc- 
tion of binary stars with intermediate circulation. Adopt- 
ing the approach of pl)| we generalize the equations for 
irrotational fluid fiow and show that stars with interme- 
diate circulation can be constructed by solving one addi- 
tional elliptic equation. We also propose a decomposition 
of the fiuid flow that represents a sequence connecting ir- 
rotational stars with corotational ones in the Newtonian 
case, and approximately corotational ones in the rela- 
tivistic case. Exploring the physical properties of these 
sequences will require future numerical work. It should 
also be possible to generalize the decomposition of the 
fiuid fiow to allow for rotation about other axes of rota- 
tion, or to modify the degree of differential rotation. Our 
approach may also be used to construct neutron stars in 
quasiequilibrium orbit about black holes. 

Relaxing the assumption of cither corotational or irro- 
tational fiow in binary stars is, at the very least, inter- 
esting as a point of principle. Such a formalism is neces- 
sary, however, for the construction of binaries containing 
rapidly rotating stars, for example millisecond pulsars, 
when the binary's orbital frequency is no longer much 
larger than the spin frequency of the star. As suggested 
by MS, approximate inspiral sequences of such binaries 
can then be constructed by fixing the circulation along 
one particular fluid trajectory, for example along the stel- 
lar equator. This can be achieved by iterating, at a given 
binary separation, over the parameter 77 in the decom- 
position of the fiuid flow until the desired value of the 
circulation has been realized. 
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Taking a divergence of the Euler equation ([37| results in 



D^H + riD'{V^K,j) = 0, 



(A5) 



APPENDIX A: AN ALTERNATIVE 
DECOMPOSITION OF THE FLUID VELOCITY 



where we have now defined 



As we have seen in Section IIVDI the decomposition 
(|52| does not necessarily lead to exactly corotational fluid 
flow in the limit = 1. An attractive alternative to this 
decomposition is 



which leads to 



-D'^+{ri-l)B' 



(Al) 



(A2) 



Retracing the derivations of Section IIV CI we find that 
equation (|56p can now be written 



H + B'Di^ 



while the continuity equation (j38[) becomes 



DH + (r/-l)A 



H + BWj^ 
Q,2 -rjBi'Bk 



B' 



(A3) 



(A4) 



= Dj{u^hBi) - D,{u*hBj). 



(A6) 



We first note that, for 77 = 0, the above equations re- 
duce to those for irrotational fiow as before. For rj = 1, 
equation (jA4[) . together with its boundary condition, is 
now solved by $ = C$, where C$ is a constant. That 
implies that, unlike our original decomposition ()52|) . the 
alternative decomposition (|Aip does reproduce corota- 
tional fiow with = in the limit 77 = 1. While this is 
a very attractive feature, this decomposition also has dis- 
advantages. Note that the term contains first deriva- 



tives of $, meaning that the term D^mj in equation (|A5|) 
will introduce third derivatives of Using the decom- 
position ([52]) . on the other hand, we avoid these third- 
derivative terms. The appearance of u^h in equation (j58p 
still introduces second derivatives of $. All of these terms 
disappear in the Newtonian limit, however, are therefore 
small, and hopefully do not affect the well-behavedness 
of the system. 



[1] T. W. Baumgarte and S. L. Shapiro, Phys. Rept. 376, 
41 (2003). 

[2] T. W. Baumgarte, G. B. Cook, M. A. Scheel, S. L. 
Shapiro, and S. A. Teukolsky, Phys. Rev. Lett. 79, 1182 

(1997) . 

[3] T. W. Baumgarte, G. B. Cook, M. A. Scheel, S. L. 
Shapiro, and S. A. Teukolsky, Phys. Rev. D 57, 7299 

(1998) . 

[4] T. W. Baumgarte, M. L. Skoge, and S. L. Shapiro, Phys. 

Rev. D 70, 064040 (2004). 
[5] C. S. Kochanek, Astrophys. J. 398, 234 (1992). 
[6] L. Bildsten and C. Cutler, Astrophys. J. 400, 175 (1992). 
[7] S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, Phys. 

Rev. D 56, 7740 (1997). 
[8] H. Asada, Phys. Rev. D 57, 7292 (1998). 
[9] S. A. Teukolsky, Astrophys. J. 504, 442 (1998). 
[10] M. Shibata, Phys. Rev. D 58, 024012 (1998). 
[11] E. Gourgoulhon (1998), gr-qc/9804054. 
[12] S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, Phys. 

Rev. Lett. 82, 892 (1999). 
[13] E. Gourgoulhon, P. Grandclement, K. Taniguchi, J.-A. 
Marck, and S. Bonazzola, Phys. Rev. D 63, 064029 
(2001). 

[14] K. Taniguchi and E. Gourgoulhon, Phys. Rev. D 66, 



104019 (2002). 

[15] K. Uryu and Y. Eriguchi, Mon. Not. R. Astron. Soc. 303, 
329 (1999). 

[16] K. Taniguchi, T. W. Baumgarte, J. A. Faber, and S. L. 

Shapiro, Phys. Rev. D 72, 044008 (2005). 
[17] K. Taniguchi, T. W. Baumgarte, J. A. Faber, and S. L. 

Shapiro, Phys. Rev. D 77, 044003 (2008). 
[18] P. Grandclement, Phys. Rev. D 74, 124002 (2006). 
[19] P. Marronetti and S. L. Shapiro, Phys. Rev. D 68, 104024 

(2003). 



8 



Erratum 

It was kindly brought to our attention by E. Goiirgoulhon that the formalism that we propose to construct binary 
neutron stars with arbitrary circulation solves only a subset of the equilibrium equations in general, but not all of 
them. 

Specifically, our formalism considers only the divergence of the Euler equation for stationary equilibrium, Eq. (13) 
in the Newtonian case and Eq. (37) in the relativistic case. However, for the Euler equation to be satisfied, both its 
divergence and its curl has to vanish. For the limiting cases of corotational or irrotational fluid flow the curl of these 
equations does vanish identically, but this is not necessarily the case for intermediate circulation. 

While our formalism still leads to valid binary initial data for fluid flow with intermediate circulation, these data 
are not necessarily in exact equilibrium for the Newtonian case or quasiequilibrium in the relativistic case. Given that 
the curl vanishes for both limiting cases, it is possible that the curl is small for intermediate conflgurations as well, 
meaning that these data may represent good approximations to equilibrium configurations. It may also be possible 
to construct equilibrium (or quasiequilibrium) data by adding another contribution to the decomposition of the fiuid 
velocity, Eq. (23) in the Newtonian case and Eq. (52) in the relativistic case, and solving for this contribution by 
forcing the curl of the Euler equation to vanish. 



